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2352134 



LOCATING MOBILE TELEPHONES 



The present application relates to the location of mobile telephone handsets. 

A number of services could benefit from the availability of accurate location 
information relating to mobile telephones. These include home zone billing, fleet 
management, emergency road services, navigation, public transport management, 
network engineering and optimisation and traffic information. 

At present, it is only possible to locate a particular handset by identifying the 
base station which is serving that handset. This will normally be the closest base 
station. In a rural area, this could cover a large geographical area. 

Existing US FCC's E91 1 requirements will, in time, require a handset to be 
locatable to an accuracy of 125m for 67% of calls and within 300m for 95% of 
calls to permit the speedy arrival of emergency services and/or breakdown services. 

Trials have shown that although accurate location can be achieved by the use 
of observed time differences, the method requires complex hardware to be added 
to the network. The method relies on the fact that GSM telephones are served by 



2 

a single base station, and measure signals from up to six neighbouring cells. It also 
relies on additional hardware to perform timing measurements of all cells in a 
network. The signals measured at the mobile from the server and neighbours will 
arrive at different times depending on distance. The absolute time of transmission 
is not known, but the time difference observed between two base stations together 
with timing measurements made on the network allow a hyperbola to be 
constructed linking points of equal time difference. With two such hyperbolae 
derived from three base stations, the intersection can be calculated to yield the 
location of the handset. Figure 1 shows this schematically. 

This approach is difficult to deploy in practice since it requires hardware 
changes to both the network and the mobiles. The changes required to the 
network involve one of two options. Either the deployment of a network of 
Location Measurement Units (LMUs) which measure the transmission time of cells 
or accurate synchronisation of the transmission times for all cells in the network. 
The latter is an unlikely choise. as it would involve costly upgrades of all base 
stations. The changes required to the mobile involve more accurate measurement 
of Observed Time Differences from surrounding cells. Furthermore, the 
mathematics for calculating the intersection is not trivial and the time differences 
will have an associated error, illustrated in Figure 2, so the hyperbolae may not 
intersect at all (Figure 2a), or may intersect at several locations (Figures 2b). As 
more time differences are employed, which might be thought to improve accuracy, 
the mathematics becomes still more complex as it is unlikely that the hyperbolae 
will intersect at a common location (Figures 2c, 2d). Figure 2e shows this situation 
most cleariy. in which the addition of hyperbola 4 has in fact worsened the 
estimated location. 

The present invention provides a method of locating a mobile telephone 
handset using the properties of the transmitted signal, such as the observed time 
differences or the received signal strength, to calculate a probability map for the 
location of the handset. This side-steps the complex mathematics involved in 



BNSDOCID: <GB 23S2134A_I_> 



3 

finding the intersections of hyperbolae. Use of the received signal strength has the 
further advantage that no hardware modifications are required to the network and 
no modifications at all are required to the mobile handset. Given this, the method 
is particularly useful in locating legacy mobiles in any situation where a positioning 
method is being deployed in a network. 

The present Invention therefore provides a method of locating a mobile 
telephone handset comprising the steps of (i) calculating a predicted signal property 
for a plurality of possible locations for the handset, (ii) comparing the predicted 
signal property with the measured signal property and (iii) deriving from the 
comparison a probability that the handset is at one or more of the said locations. 

The possible location can, for example, be distributed on a grid located at a 
likely area for the handset. That likely area can be the coverage area of the serving 
base station, or it could be an area designated as high probability by an earlier 
iteration of the same method. 

This method is highly advantageous over seeking intersections of hyperbolae 
since the computation required is simple and well suited to automation. Thus, it 
can be processed swiftly. Furthermore, it is immune to errors in the measured 
signal since only a probability is calculated as opposed to a precise location. 

Furthermore, this method can be refined significantly by combining the thus 
derived probability map with an a priori probability map, for example based on 
geographical information. This relies on the fact that it is intrinsically more likely 
that a handset will be located (e.g.) on a road or in a built-up area, as opposed to 
the middle of a desert. The combination, of the two probability maps can be a 
simple cross product of the 2D vectors, illustrated by way of example (in one 
dimension) in Figure 4. Line 20 is the varying probability derived from the method 
set out above, whereas line 22 is a step function inferred from local geographical 
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data, in this case a road passing through a barren area, 
position of the handset is at 24. 



Thus, the most likely 



An alternative a priori probability map is one based on historical data, 
possibly combined with geographical data. Thus, a peak probability which moves 
along a path parallel to a road indicates that the handset Is on that road. This is 
particularly so if the peak probability is moving at speed. A bend in the road would 
also give information as to the location along the road. Alternatively, a long rest 
period adjacent a spot feature such as a congested junction offers a refinement to 
an instantaneous position. Previous position data also enables a corrective 
displacement relative to the peak probability to be calculated. 

The preferred signal properties are time differences and signal strength, 
particularly signal strength for the reasons given above. Where signal strength is 
employed, it is also preferred that pre-calculated signal strength arrays are loaded 
for the serving cell and any measured cells involved. 

Several probability maps could be combined simultaneously or sequentially. 
For example, the maps based on time differences, geographical data, historical data 
and signal strength could all be crossed to identify a peak probability. 

An embodiment of the present invention will now be described by way of 
example, with reference to the accompanying drawings, in which:- 

Figure 1 is a schematic illustration of the intersecting hyperbolae method; 

Figure 2 is a schematic illustration of errors in the intersecting hyperbolae 
method; 

Figure 3 illustrates an example of a probability map; 
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Figure 4 shows the combination of probability maps; 

Figure 5 shows a part of a GSM network; and 

Figure 6 illustrates time differences in signal arrival times. 

B^clcground to GSM networks 

In a GSM network mobiles are only served by one base station but they 
always measure signals from neighbouring cells (typically six neighbours) (see 
Figure 5). This signals from these base stations all arrive at different times 
depending on when they were transmitted and how far away they are from the 
mobile. This document will conform to the following naming convention for these 
quantities. They are: 

OTD: Observed Time Difference. This Is the time difference between the 
arrival of the signal at the mobile from the serving base station and the arrival of 
the signal from a neighbouring base station. This is illustrated in Figure 6. 

RTD: Real Time Difference. This is the difference between the transmission 
of a frame from one base station (BS) and the transmission of a frame from 
another. If BS1 transmits a frame at time t = 5 and BS2 transmits a frame at time 
t = 7 then the RTD between them is 7-5 = 2. If all the base stations In a network 
transmit their frames at the same time then the network is said to be synchronised 
and all the RTDs are zero. 

GDT: Geometric Time Difference. This is the part of the time difference of 
the arrival of two frames due to geometric separation of the sites and the mobile. 
This difference in time taken for a transmission to travel from BS1 to the mobile 
and the time taken for a transmission to travel from BS2 to the mobile. 
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The OTD for any pair of sites is the sum of the RTD and the GTD 
(OTD = RTD + GTD). 

The mobile measures the OTD for the serving base station and each of the 
neighbours. The GTD is what is required to calculate the position of the mobile. 
This can either be obtained by subtracting the RTD from the OTD or by using a 
synchronised network such that the OTD and the RTD are the same. In this 
description, it is assumed that the networic is synchronised, but the correction is 
straightforward. It is also assumed that propagation from the base stations is line 
of sight. This assumption allows us to convert the time difference to a distance 
difference by dividing by the speed of light. 

Location method 

As set out above, location of mobile handsets according to preferred 
embodiments of the present invention can be by use of observed time differences 
or signal strength measurements. In this example, the detailed explanation of the 
method will be given mainly in relation to observed time differences. However, the 
principal as applied to measured signal strengths is identical and the resulting 
mathematics equally applicable. 

Instead of constructing a hyperbola for each OTD measurement, a preliminary 
estimate of the location of the mobile is made (e.g. using the cell of origin) and then 
a grid Is constructed around the estimated position. Then for each square of the 
grid, a prediction is made of what the OTD or signal strength measurements would 
be, for all the measured base stations, if the mobile were in that square. In 
practice, the signal strength for relevant areas around each site would be pre- 
calculated so that the process would involve loading precalculated arrays of 
predicted, signal strength instead of a more time-consuming signal strength 
calculation. 
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Those predictions are compared with the measured values to give a figure 
of probability of the mobile being in that square. A graphical representation of a 
probability map is shown in Figure 3. The darker the shading, the higher the 
probability of the mobile being in that square. The resolution of the grid could 
either by very fine - say 1 m or 5m, or to reduce the calculations required, the 
process could be done in stages with an initial coarse grid {say 1 00m) followed by 
smaller and finer grids centred on the high probability area. 

Further advantages can be added by combining the information from the 
probability map with information about the area concerned. This might include 
high-resolution digital mapping data information or corrections for areas where 
propagation is known to be non-line-of-sight. Alternatively, the a priori map could 
simply be an averaged historical map of handset location. 

Since some of the calculations involved in calculated the predicted OTDs, 
which are compared with the real measurements, do not change over time, it may 
be possible to speed up the location calculation by doing a certain amount of 
calculation in advance and storing it for later use. 

The remainder of this description outlines the numerical method for 
computing the location of a handset based on the measured delay between signals 
received from several sites. The only restriction on the method is that the 
measured delays are statistically independent. A simple example is the situation 
where all delays are relative to the signal received from one principal site. The 
method is also applicable to the use of received signal strength measurements 
(rxLev) made on the handset. These measurements are assumed to be available 
through an interface to the BSC (Base Station Controller). The method would be 
the same if a mobile were providing the information via a short message. 

In order to obtain an estimate for the handset location, one must also make 
assumptions regarding the errors in the measured delays or signal strengths. The 
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following assumes that these errors are normally distributed, but further 
experimentation may show that an alternative distribution should be adopted. 

In the following; 

lower case letters are used to signify continuous random variables, e.g. x; 

a lower case p{) signifies a probability density function, e.g. p(x); 

upper case letters are used to denote a tiny range of values of a continuous 

random variable, e.g. 

X refers to values of x lying in the range [X-6X, X + 6X1 
Y refers to values of y lying in the range IY-QY, Y + 6Y1; 
an upper case PO signifies a probability. This is invariably obtained by 

integrating a probability density function between appropriate limits: 



x*ex 



x-sx 



We shall write integrals of the above form using the more concise notation 

P{X)= \p{x)dx 

X 

For simplicity of presentation it will be assumed that there are just two 
independent time difference (OTD) measurements (T,.!^) or corresponding signal 
strength measurements. The generalisation to higher numbers of independent 
measurements is trivial. Note that the measurements have been written in upper 
case. This signifies that they refer to finite ranges [T,-5T,.T, + 5t,] and ITj- 
5T2.T2 + 5Tj1 of continuous random variables t, and tz- The widths (26T, and 
25T2) of these ranges are determined by the accuracy of the measuring equipment. 



It is necessary to calculate: 

/»(A'.r|r„r,) 



{1) 
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This is the probability that the handset Is at {X,Y) given that we have timing 
measurements (Ti.Tj) or corresponding signal strength measurements. 



We can use Bayes* Theorem (i.e. that P(A,B> = P(A|B) = P(B|A)P(A) to 
write (1) as 

The terms on the right-hand side can be written as integrals of the joint 
probability density functions p(x,y,t,,t2). 

J J J j p{x,y,t„t^)dxdydt,dt^ 
P{X,Y\T,,T,)= '' : 1 (3) 

Tf Tj -«o -co 

We can use Bayes' Theorem again to write P(x,y,t,,t2) in terms of things that 
can be calculated: 

J J J j P{f„t2\x,y)p{x,y)dxdydt,dt^ 
j J J j p{t„t2\x,y)pix,y)dxdydt,dt^ 

7, Tj -«> 



The factor p(x,y) is simply the probability density function for a spatial 
location. The simplest assumption that could be made is that p(x,y) is uniformly 
distributed over the area we are considering. Other choices are possible and may 
give improved results, for example, areas of high traffic could be given a higher 
value than areas of low traffic. 
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Assuming that p(x.y) is uniformly distributed over the mobile's serving cell, 
(4) can be simplified by cancelling the p(x,y) terms, and replacing the spatial 
integral in the denominator by an integral over the serving cell. 



p{x. r 1 7-,, r,) = yj'^p(/../,ix.>.WMJr, 



(5) 



T, r, Serymg 
C«0 



The factor pd^t^jx.y) is the probability density function for the timing 
measurements (or corresponding signal strength measurements) given that the 
handset is at (x,y>. This can be calculated if some assumptions are made about the 
errors in the timing or signal strength measurements. The locations of the 
transmitting sites are known, and for any mobile location (x,y) the ideal time 
measurements (or corresponding signal strength measurements) can be calculated 
that would arise in the absence of error. These can be called "ideal measurements" 
t;(x.y) and t;(x.y). Assuming also that the errors in the measured times are 
normally distributed with zero mean (but possibly different variances O, and Oj). 

Since the measurements are statistically independent, we can write 



TTtT, 



-exp 



2<Ti 



(6) 



Now if our timing or strength measurements are made with sufficiently high 
accuracy (5t, < <0„ Stj < <02) the integrals dt,dt2 can be approximated in (5) to 
give 

1 .-r.Ml_i_.xp(-^ M.* 



p(x.y\T,.T,)': 



CM 



(7) 
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At a small enough spatial resolution, the remaining integrals in (7) can be 
approximated in a similar manner. This Is essentially what the following algorithm 
does. 



Algorithm for determining handset position 

a) Scan through an area with fixed step sizes in the x and y directions. These 
should be sufficiently small. 

(b) At each pixel {X,Y), calculate the "ideal measurements" t;(X,Y) and t2(X,Y), 
and then evaluate 



-Jin a. 



exp 



exp 



r Ml 



2<Ti 



<8) 



(c) The probability that the handset is in a pixel (X,Y) is simply 

P(YY\TT^'^ (9) 



Fixeis in 
Servtng Cell 



Note that (9) is an approximation to (7) but is valid if the pixels are 
sufficiently small. Clearly, the smaller the pixel size, the better the 
approximation. 

(d) Once P(X,Y|T,T2) has been worked out at every pixel in the search area, the 
most likely location for the handset is given by the expectation 

ServtngCeU ( 1 Q) 

e) The variance (if calculated) will give an idea of the accuracy for the 
prediction (i.e. how well localised it is). 
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JX-E{xy\ (X-E{xy\ 



Ser\'mg Cell 



(11) 



This could be used to draw a bounding circle around the point 



Alternatively, if all calculations are stopped at step (b> and the largest V{X.Y) 
in the search area issued as an indication of the handset position, then the same 
result could be obtained by sinnply taking the mininnunr> of 

\2 /. 



fzkziY 
^ ) 



I.e. use X' as a maximum likelihood estimator for (X,Y). 

On an Implementation level it is assumed that the relevant pre-calculated 
signal strength arrays could be downloaded to the BSC or other mobile location 
platform on a regular basis. 

The present invemion therefore offers a numerical method for calculating the 
location of a handset based on independent timing or signal strength 
measurements. The method attempts to determine the most likely position for a 
handset and is able to give an idea of the accuracy of the result. The method has 
been described above in a general form and the final implementation would allow 
the user to "tune" the method by varying the following parameters: 

Pixel size. This needs to be sufficiently small. 

Distribution functions for the time delays between pairs of sites or signal 
strengths from different cells. It is possible to experiment with different 
standard deviations. 

p(x.y) . This acts as a weighting function when calculating the V values. For 
example, the calculation could be biased to favour areas of high traffic. 
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The method proposed herein is computationally expensive compared to other 
techniques of handset positioning (based on intersecting hyperbolae for example), 
but is correspondingly more straightforward to automate. It also makes use of only 
one set of timing measurements (or corresponding signal strength measurements) 
and so it does not benefit from previous measurements/calculations that may be 
available to the system. Prior knowledge of this sort can be used to obtain faster 
(and more reliable) estimates of handset position. 

It will of course be apparent to the skilled reader that many variations could 
be made to this embodiment without departing from the scope of the present 
invention. 
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CLAIMS 

1 A method of locating a mobile telephone handset comprising the steps of ; 

(i) calculating a predicted signal property for a plurality of possible 
locations for the handset, 

(ii) comparing the predicted signal property with the measured 
signal property and 

(ill) deriving from the comparison a probability that the handset is 
at one or more of the said locations. 



2. 



A method according to claim 1 in which the thus derived probability map is 
combined with an a priori probability map. 

3. A method according to claim 2 in which the a priori probability map includes 
geographical information. 

4. A method according to claim 2 or claim 3 in which the a priori probability 
map includes historical data. 

5. A method according to any one of claims 2 to 4 in which the combination 
of the two probability maps is a product of the two probabilities at each 
location. 

6. A method according to any one of claims 2 to 4 in which there are several 
probability maps each based on a different factor. 

7. A method according to claim 6 in which the maps are combined sequentially. 

8. A method according to claim 6 in which the maps are combined 
simultaneously. 
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9. A method according to any preceding claim in which the signal property is 
the measured signal strength. 

10. A method according to any one of claims 1 to 8 in which the signal property 
is an observed time difference. 

11. A method according to any one of claims 1 to 8 in which a probability map 
is also generated based upon observed time differences. 

12. A method of locating mobile telephones substantially as described herein 
with reference to and/or as illustrated in the accompanying drawings. 

1 3. Apparatus for locating mobile telephones adapted to implement a method as 
set out in any one of the preceding claims. 
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